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Abstract 

Thin- layer Navier- Stokes calculations 
for wing-fuselage configurations from sub- 
sonic to hypersonic flow regimes are now 
possible. However, efficient, accurate 
solutions for using these codes for two- 
and three dimensional high-lift systems 
have yet to be realized. A brief overview 
of salient experimental and computational 
research is presented. An assessment of 
the state-of-the-art relative to high-lift 
system analysis and identification of is- 
sues related to grid generation and flow 
physics which are crucial for computation- 
al success in this area are also provided. 
Research in support of the high-lift ele- 
ments of NASA's High Speed Research and 
Advanced Subsonic Transport Programs which 
addresses some of the computational issues 
is presented. Finally, fruitful areas of 
concentrated research are identified to 
accelerate overall progress for high lift 
system analysis and design. 


I . Introduction 

An area of special interest to aerospace 
designers is high-lift systems. Future 
transport aircraft will have multiple re- 
quirements playing important roles in 
their design. These requirements include 
improved energy efficiency, reduced noise, 
and lower maintenance costs. Improved 
high-lift concepts for subsonic transports 
may result in designs which have increased 
section thicknesses, larger aspect ratios, 
lower sweeps, optimized multi-component 
designs, highly integrated propulsion sys- 
tems, and integrated pneumatic concepts 
such as circulation control. Conversely, 
transports designed for supersonic cruise 
typically have geometric characteristics 
(highly swept, slender wings) which do not 
lend themselves to efficient aerodynamics 
at low subsonic speeds and moderate-to- 
high angles of attack (flight conditions 
associated with takeoff and climb-out) . 
The need for high-lift augmentation con- 
cepts is further accentuated by contem- 
porary community noise standards and traf- 
fic congestion at Air Traffic Control sta- 
tions. While there is ongoing research 
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in this area at national research laborat- physics and geometry, 
ories and private industry, (Brune and 

McMasters 1 provide an extensive review of High-lift systems for supersonic config- 
computational high-lift design in practice urations differ from those for subsonic 

at industry) , a critical need exists for transport systems both in system geometry 

further innovations. For this to be real- and physics which dominate the flow. In 

ized, major breakthroughs in several areas order to achieve the desired high levels 

must occur. From a computational perspec- of supersonic cruise efficiency, many ad- 

tive improved methods are needed to analy- vanced supersonic configurations employ a 

ze geometrically complex systems and in- low aspect ratio, highly swept wing, Un- i 

elude key physics, such as flow separa- fortunately, these configurations typical- 

tion, transition and turbulence, which ly have poor low-speed performance charac— 

dominate the flow fields. This paper at- teristics. Often, the low speed perfor- = 

tempts to review some of the issues which mance characteristics of these systems are I 

are crucial for computational fluid dynam- enhanced either by attached flow or vortex 

ics (CFD) to truly complement ground and flaps along the wing leading edge. At- 

flight based research and development for tached flow flaps are designed to suppress j 

advanced high-lift systems. the formation of leading-edge vortices. : 

Conversely, vortex flaps are designed to 

Advanced transport designs currently position the leading edge vortex within E 

receiving considerable attention include the bounds of the flap chord to provide a = 

configurations designed for supersonic component of thrust which results from a 

cruise, such as the High Speed Civil vortex induced suction force. The trail- = 

Transport (HSCT) , as well as more conven- ing edge flap system for these configura- 
tional subsonic transports. High-lift tions may consist of a segmented system of 

systems for subsonic transports, typically hinged flaps. Figure 2 shows the 

use deflected leading edge slat surfaces schematic of a low aspect ratio highly 

and trailing edge slotted flaps for lift swept wing configuration tested at NASA ; = 

augmentation, see Figure 1. Figure la the leading edge flap segments can be 

shows a relatively simple 3-component sys- deflected independently about the hinge 

t em consisting of a main element, slat and line. The trailing edge flap segments can 

single slotted flap configuration which also be deflected independently about the = 

was tested by Lockheed-Georgia 2 . A more flap hinge lines. Also shown in Figure 2 

complex system shown in Figure lb depicts are schematics of attached flow and vortex 

a double slotted trailing-edge flap in flap concepts. Grid systems to model 

addition to the slat and main element. these complex, segmented geometries must 

This configuration was tested in the NASA be highly versatile. In addition, the 

Langley Low Turbulence Pressure Tunnel computational methods employed in the 

(LTPT) . Subsonic high-lift systems, when study of high-lift systems for supersonic |§ 

fully deployed, can have regions of configurations must be capable of captur- 

separation on the slat or flap, or in the ing vortex structures with minimum smear- ^ 

cove regions. The subsonic systems may ing and phase distortion since the nature 

also have confluent boundary layers as a of the flow is highly vortical in these 

result of the strong interaction of shear systems, 

layers. The subsonic high-lift systems 

are thus viscous dominated flow fields and The problems of high-lift system analys- 

have considerable geometrical complexity. is often include such issues as engine „ 

Computational methods for high-lift sys- airframe integration and three-dimensional 

terns must be carefully chosen to incor- effects resulting from wing sweep, pylons, — 

porate these varying requirements in flow partial— span— flap deflections and tip ef— _ , 



II .Literature Overview 


fects. These attributes yield extremely 
complex geometries and attendant complex, 
interactional physical phenomena. While 
incremental progress by way of interac- 
tional methods or in development of quasi- 
three-dimensional analyses is being accom- 
plished, we are not confident that such 
progress will computationally support in- 
novative breakthroughs for future high- 
lift system design and development. Low 
cost, computationally efficient solutions 
(multi-grid, local time stepping) for 
three-dimensional steady flows for rela- 
tively simple geometries are available 
using thin-layer Navier-Stokes codes. We 
feel emphasis on these methods will yield 
substantial progress toward alleviation of 
the two major obstacles to accurate, effi- 
cient high-lift system analysis, viz., 
complex geometry and physics. Hence, our 
principal perspective will focus on Navi- 
er-Stokes solutions . In this paper, we 
explore some of the crucial issues which 
must be successfully addressed in order to 
develop computational methodology to 
analyze three-dimensional high-lift sys- 
tems . 

The following sections present a brief 
overview of experimental efforts useful 
for code calibration or validation along 
with a brief assessment of the computa- 
tional state-of-the art. Geometric con- 
siderations are addressed within the 
context of the implications for grid 
generation. Issues related to code 
algorithms and dominant physics are 
briefly discussed for high-lift system 
applicability. Some of the salient 
efforts currently being pursued at NASA 
Langley are also presented. Concluding 
remarks consist of suggestions of major 
research areas where coordinated work is 
required to sustain progress for 
computational high-lift system analysis 
and design. 


Experimental database 

While it is not the intent of this paper 
to review the available high-lift research 
literature, a brief overview of some sali- 
ent reports is appropriate. Among the 
published data on multi-element airfoils, 
Braden, Whipkey, Jones and Lilley 2 report 
on a study of the confluent boundary layer 
development and separation characteristics 
on a NASA GAW-1. The section was equipped 
with a 29% chord single-slotted trailing 
edge flap and a 15% chord leading edge 
slat. Various combinations of slat and/or 
flap deflections and angle of attack were 
investigated in the study. The report con- 
tains surface pressure measurements on the 
airfoil as well as lift versus angle of 
attack curves for various flap/slat 
arrangements. Surface oil flows were used 
to provide flow visualization of boundary 
layer transition patterns* Boundary layer 
velocity profiles, turbulence intensities 
and Reynolds shear stresses for the con- 
figurations are reported under a separate 
cover . A supplement to this report con- 
tains over 30000 sets of laser velocimetry 
(LV) derived boundary layer and wake data 
for the various combinations of geometric 
arrangements and angles of attack. In 
addition, off body flow field data were 
obtained using hot wire and LV. 

Valerazo, Dominik, McGhee, Goodman and 
Paschal 3 have conducted multi-element air- 
foil optimization studies for maximum 
lift. This is a cooperative study between 
NASA and Douglas Aircraft Company. The 
primary focus of the study was to discern 
the high Reynolds number sensitivities of 
the multi-element airfoils at chord 
Reynolds numbers up to 16 million. The 
high-lift system consists of a double 
slotted flap and a single slotted slat as 
shown in Figure 2b. Among the data that 
are presented is the variation of C, 
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with Reynolds number, the variation of 
c imax with slat and the effect of flap 



gap on C lmax . No flow field surveys are 
reported from the study. In addition, 
there were no mechanisms used in the study 
for transition detection. 

Olson and Orloff 6 report on the study 
of an airfoil with flap arrangement con- 
ducted in the NASA Ames Research Center 
7 by 10 foot tunnel for Mach Number of 
0.06 and a Reynolds number of 1.3 million. 
Surface pressure measurements, Reynolds 
stresses and detailed measurements^of mean 
velocity in the boundary layers, wakes and 
merging layers are reported. The data 
should be considered purely incompressible 
and codes with compressible formulations 
will have difficulty in simulating this 
extremely low Mach number. 

Wentz, Seetharam and Fiscko 7 tested an 
aileron and a fowler flap applied to a 
GAW-1 airfoil. The experiment was conduc- 
ted at M = 0.13 and Reynolds number of 2.2 
million. Aileron control effectiveness 
and hinge moments are presented for vario- 
us gaps from 0 to 2% chord. For the 
fowler flap study, pressure distributions 
for various flap settings were obtained 
for a limited angle of attack range. 

Adair and Horne 8 present pressure and 
velocity characteristics in the vicinity 
of the flap of a single slotted airfoil at 
a Mach number of 0.0 9 and a chord Reynolds 
Number of 1.8 million at the NASA Ames 
7x10 foot tunnel. They report strong con- 
fluence effects on the boundary layer dev- 
elopment on the flap suction surface due 
to the presence of a strong jet emanating 
from the slot flow. The flap is separation 
free except at the trailing edge where 
intermittent separation is observed. As 
a result of this the data may only be j?f 
limited use for steady state calculations. 
The flap wake development is reported to 
be asymmetric due to strong destabilizing 
curvature effects on the suction side. 
These data should therefore provide some 
guidance for studies of non-equilibrium 
effects on turbulence . 


9 10 

Morgan and Morgan and Paulson report 

on the study of static longitudinal and 
lateral directional aerodynamic 
characteristics of an advanced aspect 
ratio 10 and aspect ratio 12 supercritical 
wing transport model. The model was 
equipped with a high-lift system 
consisting of a full-span leading edge 
slat and partial-span and full-span 
trailing edge flaps. The Reynolds number 
of the tests varied from 0.97 to 1.63 
million over a Mach number range of 0.12 
to 0.20. The model was tested at angles 
of attack from -4 to 24 degrees and 
sideslip from -10 to 5 degrees. The model 
has engine nacelles, landing gear and 
movable horizontal tails. Six basic wing 
configurations were tested. These 
consisted of cruise (nested case) , 
partial-span flap, full-span flap, full- 
span flap with low-speed ailerons and 
full-span flap with high-speed ailerons 
with slat and flap deflected to represent 
takeoff and landing conditions. Lift, 
drag and pitching moment data are present- 
ed for various cases. 

Nakayama, Kreplin and Morgan 11 report 
detailed flow field measurements for a 
three-element airfoil with a conventional 
slat and single slotted flap. Reynolds 
stress distributions an 4 mean flow 
measurements are presented on the main 
element and in the flap and wake regions. 
These suggest strong confluence effects in 
the flap region involving a jet-like 
stream from the flap-airfoil gap, the wake 
region of the main element with slat and 
the boundary layer on the flap itself. 

The above cited works provides some data 
for code calibration or validation. 
However, there are large voids in the data 
base and measurements are not in suffi- 
cient detail to understand the complex 
flow physics that a computational study 
seeks to model. The flow in the gap 
region and pressure sides of many of these 
configurations needs to be documented 


fully. There is also a need to obtain 

* f * t 

turbulent fluctuations data (u , v , u v ) , 
transition location and mapping of flow 
confluence. Further studies in this area 
must be designed to closely follow the CFD 
needs to construct proper turbulence 
models and flow modules for Reynolds 
averaged Navier-Stokes calculations . 

Computational database 

Brune and McMasters 1 provide an excel- 
lent review of existing computational met- 
hods for analysis of high-lift systems. 
The status of these methods can be summar- 
ized as follows: there are presently no 
truly three-dimensional CFD methods for 
high-lift studies. Most three-dimensional 
studies use quasi-three-dimensional vis- 
cous approaches such as three-dimensional 
inviscid codes coupled with two-dimension- 
al boundary layer codes. Existing vis- 
cous, two-dimensional airfoil codes can 
be classified, according to Ref .1, as, (I) 
Coupled Attached-Flow Methods, (II) Coupl- 
ed Separated-Flow Methods, (III) Navier- 
Stokes Methods and (IV) Design and Optimi- 
zation Methods. Categories (I), (II) and 
(IV) are widely used in industry today, 
while Category (III) is considered to be 
at the developmental stage. Both the at- 
tached and separated flow methods (Catego- 
ry (I) and (II) ) are based on interaction- 
al boundary layer approaches while the 
Design and Optimization methods are clear- 
ly a patchwork of methods (I), (II) and 
simple inviscid analyses. In category 
(I), a boundary layer method is coupled to 
an inviscid flow calculation 12,13,14 and 
in (II) , some form of modelling of the 
separated region is attempted 15 ' 1 ' 17 . The 
attached flow methods provide good agree- 
ment for lift at low angles of attack, 
where there is no flow separation. The 
separated flow models have been successful 
for some cases to compute maximum lift up 
to stall. These methods are at best use- 
ful in a limited fashion and do not 
promise to provide a successful methodol- 
ogy for high-lift system design. There 


are also a few applications of two-dimen- 
sional Navier-Stokes solvers for high-lift 
configuration analysis in the literature. 
Schuster and Birckelbaw 18 and Shima 19 have 
obtained two-dimensional Navier-stokes 
solutions for multi-element airfoil sys- 
tems using patched structured-grids. 
Using an unstructured-grid solver, Mav- 

i i 20 

riplis and Martinelli have also obtained 
solutions of two-dimensional multi-element 
airfoils. This work may well be a 
bellwether for high-lift system 
computations . 

Ill . Geometrical Considerations 

Complex geometry issues associated with 
high-lift system analysis are non-trivial 
to say the least. Even when the problem 
is simplified to a wing with deflected 
surfaces (disregarding pylons, engines, 
flap track fairings, etc.), the task of 
surface modelling and field discretization 
is formidable. Geometries which are dis- 
continuous in the streamwise and spanwise 
directions offer a significant challenge 
to the CFD community. Within these dis- 
continuous regions flow interactions are 
occurring which can have significant and 
dominant effects on the resulting flow 
field. An example of this is the vortical 
flow occurring at the edge of a partial 
span leading-edge flap as it is deflected 
on a highly swept wing. The following 
sections address in some detail the manner 
in which CFD code developers are address- 
ing these issues. From the structured 
grid perspective, single block, multi- 
block and Chimera schemes are each addres- 
sed. The promising work going on in the 
development of efficient unstructured grid 
generation techniques is also discussed. 
Finally, zonal methods are addressed in- 
cluding an example of their applications. 

Structured-grid solvers 

(a) Single and Multi-block methods 


The rapid progress in CFD of the last 



decade has made it possible to analyze 
simple wing-body geometries with relative- 
ly little effort. This is due, in part, 
to efficient grid generation techniques 
and acceleration techniques such as mesh 
sequencing, local time stepping and multi- 
grid techniques. Structured-grid al- 
gorithms, such as TLNS3D 21 ' 22 ' 23 and 
CFL3D 24 , have shown that for many steady 
flow problems, efficient solutions are 
possible using multigrid acceleration 
schemes. However, for the multi-element 
problems, the single block structured sol- 
vers are difficult to use. Fortunately, 
multi-block versions of these codes are 
currently being developed. These multi- 
block solvers may have the power to 
analyze complex domain problems by break- 
ing the flow domain into smaller sub- 
domains or blocks of individual grid 
topology {such as grid system for each 
component of a multi-component system) . 
The appropriate set of flow equations in 
each of these blocks can then be solved. 
Another significant development is the 
availability of powerful new grid genera- 
tion packages which in the hands of 
experienced users can be used to do vir- 
tually any type of gridding (C-0; C-H; C- 
C) with relative ease. Among the most 
promising grid-generation packages are 
GRXDGEN 25 and EAGLE 26 , These are both 
user-friendly packages for generating two- 
dimensional and three-dimensional struc- 
tured volume grids for finite volume 
analyses. Single or multi-blocked grids 
may be generated using these packages. 
The grid systems that may be constructed 
in the multi-blocks may or ma^ not have C 
(common grid points) or C continuity 
(slope continuity as well as common grid 
locations) at the interface of these 
blocks. Depending on the nature of these 
interface conditions many variations of 
boundary coupling between various blocks 
are possible. 

Previously, the utility of using 
structured-grid solvers for multi-element 
airfoil cases has been explored in the 


context of two-dimensional flows. Using 
a st ructured-qrid solver as the base code, 
Schuster and Birckelbaw developed solu- 
tions for the multiple element airfoil 
problem by a multi-block approach using 
two-dimensional Navier-Stokes solutions . 
Figure 3 shows the schematic of the multi- 
ple-block grid topology used by them for 
a two-element airfoil. Figure 3a shows 
the arrangement of the various blocks in 
physical space and Figure 3b shows the 
arrangement in computational space. The 
line marked S is a line of singularity 
where all three blocks intersect and it 
requires special connectivity relations. 
The composite grid in Figure 3c is ob- 
tained by an iterative approach such that 
the grid systems in the regions retain C 
continuity at the block interfaces. The 
flow solver used in Ref, 18 is a modified 
ADI scheme closely related to the Beam- 
Warming algorithm and the turbulence model 
used in the calculations is the Baldwin- 
Lomax model. Schuster and Brickelbaw ob- 
tain a reasonable comparison of C x with 
experimental data at angles of attack up 
to stall as can be seen from the lift ver- 
sus angle of attack curve shown in Figure 
3d. They also state that their solution at 
the stall angle and beyond did not con- 
verge to a steady state solution. The 
curve shown by the dotted line in the fig- 
ure is an average of the oscillatory solu- 
tion. The C p predictions on the main ele- 
ment obtained by them (not shown here) 
indicate some systematic variations from 
experimental measurements, the cause of 
which was unknown. 

Shima 19 also obtained Navier-stokes 
solutions for a multi-element airfoil sys- 
tem using a patched grid system. The grid 
generation for the multiply connected 
domain in this work is again a non-trivial 
problem. Here, the composite grid is 
obtained in a two-step process. Initial- 
ly, a potential flow solution around the 
multi-element airfoil is generated using 
a panel method. Next, conventional grid 
generation techniques 2 , using finite-dif- 


ference methods, are employed where the 
computational co-ordinates are now the 
known potential and st reamf unctions around 
the multi-element airfoil. This allows 
control of grid spacing required for the 
Navier-Stokes solutions near the body. 
The flow solver in Ref. 19 uses an upwind 
(Total Variation Diminishing or TVD) sche- 
me modified for low Mach number applica- 
tions. Computed solutions are compared to 
experimental measurements of Foster, Irwin 
and Williams . The results obtained for 
a two-element configuration (consisting of 
main element and flap) are shown in Figure 
4 reproduced from Reference 19. The agre- 
ement between experimental data and com- 
putations for C 1 is reasonable; the stall 
angle of attack predicted from the solu- 
tion by the averaging method similar to 
that used in Reference 18 is under-predic- 
ted in the calculations. The authors pos- 
tulate that this could be a result of num- 
erical problems. The results for stall 
and post-stall cases are once again suspe- 
ct since they are obtained by averaging an 
oscillatory solution obtained by the com- 
puter simulation. 

While these results for two-dimensional 
cases suggest the utility of multi-block 
systems for high-lift analysis, further 
research is required to establish the usa- 
bility of such methods for a highly com- 
plex three-dimensional configuration. For 
three-dimensional applications, the multi- 
block methods with rule based expert sys- 
tems may provide a natural way for genera- 
ting structured grids for analysis. Dan- 
29 

nenhoffer discusses the development of 
such a system for two-dimensional multi- 
body configurations. With elements in 
close proximity, the nature and quality 
of such grids and their resultant sensiti- 
vity to overall flow solution need to be 
examined closely. Many of the finite vol- 
ume structured solvers are highly sensi- 
tive to grid quality. The inability of 
these methods to provide reasonable simul- 
ations in regions where the grid may be 
highly distorted and stretched, such as. 


in the narrow regions of flow passage aro- 
und the multi-element airfoil case, may 
limit the use of these methods for 3D- 
high-lift analysis . The complexity of the 
grid-generation and flow solver may also 
have some bearing on their eventual accep- 
tance . 

(b) Overlapping Grids/Chimera Schemes 

In addition to the multi-block method, 
grid overlapping methods are another com- 
monly used technique for domain decomposi- 
tion. In overlapping schemes the sub- 
domains and the grid systems associated 
with them may overlap, or it may be pos- 
sible to embed one sub-domain completely 
in another. In the "chimera scheme 1 ' 30 , 
the regions of a grid common to others is 
removed thereby creating voids or holes 
inside the grid. Baysal et. al. , have 
looked at the quality of chimera solutions 
by studying the solutions with and without 
embedding for a test problem and conclude 
that there are only "minor" differences 
between the solutions. If this is true, 
chimera schemes may offer the flexibility 
to study multi-element airfoil flows. An 
example of a chimera grid developed at 
NASA Langley Subsonic Aerodynamics Branch 
for the GAW-1 airfoil with a deployed slat 
is shown in Figure 5 . Figure 5a is an 
example of a sub-domain which consists of 
the slat geometry. Each sub-domain (see. 
Figure 5a, b) contains a "hole" or void in 
it which is a region of overlap of another 
sub-domain. The void is identified for 
each sub-domain in a preprocessing step. 
The solution strategy for the composite 
flow field involves computation of flow 
fields in each sub-domain with the associ- 
ated boundary values including those for 
the boundaries of the void region. Since 
the boundary values for the voids are gen- 
erated iteratively (by solutions from the 
sub-domains that create the voids) , 
convergence of these methods depends 
strongly on how well the boundary values 
are approximated. 



Figure 6 shows Euler solutions obtained 
by Biedron 32 using CFL3D employing an 
overlapped grid option for an airfoil with 
a slat. The calculated conditions are at 
M=0 . 5 and a - 7.5°. These excellent 

results suggest that the overlapped grid 
option may be exploited to generate grid 
structures over multi -component airfoils. 

An important advantage of the chimera 
scheme or the overlapped grid methods is 
the relative ease with which structured-^ 
grids can be generated around "simple" 

sub-domains of a complex three-dimensional 
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domain. Buning, Parks, Chan and Renze 
describe the application of a chimera 
scheme for the space shuttle ascent 
geometry. The component grids were 
generated using a hyperbolic grid 
generation technique which is faster than 
elliptical grid generators. Due to the 
complexity of the geometry, the grid 
joining process does become somewhat 
involved at the intersection of 
geometrical components. Further 

innovation in the form of "collar grids" 
were required to develop solutions for the 
shuttle ascent geometry. An example of a 
"collar grid" for a cylinder intersecting 
a curved surface is shown in Figure 7 (re- 
produced from Ref .34) . Figure 7a shows the 
combined collar surface grid. The white 
region in the figure is the void in the 
cylinder and the plane surfaces. The col- 
lar grid separates the intersecting sur- 
faces and acts as a transitional zone bet- 
ween them. Figure 7b shows a slice of the 
collar grid and the chimera grids around 
it* 

The overlapping schemes and in 
particular, the chimera scheme provide a 
simple way to generate computational 
grids. However, further study is required 
to sort _out any sources of error in such 
an approach before recommending these 
methods as a panacea for high-lift system 
analyses. Buning et.al . 33 point out that 
while the accuracy of their solutions im- 
proved with improved modelling of the geo- 


metry, the accuracy required for wing 
loading analyses is significantly higher 
than obtainable with chimera schemes . 
Effective use of chimera is also limited, 
according to them, by difficulties for the 
scheme in implementing turbulence models 
based on length scales for multi-body con- 
figurations . 

Unstructured-arid methods 

Navier-Stokes solvers using unstruc- 
tured-grids (triangular, tetrahedral mesh- 
es) are relative newcomers to the field. 
While finite element methods using trian- 
gular and quadrilateral and tetrahedral 
elements have been used in the past, their 
applications have been limited to low Rey- 
nolds number flows. Mavriplis , Mav- 

2 0 

riplis and Martinelli and Mavriplis and 
Jameson 36 have led the way in developing 
viable solutions to flow over airfoils. 
The implementation of multigrid strategy 
and turbulence modelling for the two-dime- 
nsional cases are major assets to the flow 
solver used in Reference 20. At present 
several multi “Component flows have been 
calculated using this version of the code 
and good agreement with experimental data 
has been obtained for many of these cases. 
For example, Reference 20 documents 
solutions of multi-element airfoils which 
show excellent agreement with experimental 
data using a two-equation (k-e) model. 

The advantage of using unstructured- 
grids for analyzing high-lift systems is 
obvious. They are capable of properly 
modelling all the geometric complexity as- 
sociated with high-lift systems in a 
straightforward manner . Figure 8 shows an 
unstructured-grid system developed at Sub- 
sonic Aerodynamics Branch (SAB) , NASA Lan- 
gley Research Center, for the study of a 
fully deployed low speed multi-element 
airfoil. An important merit of these 
methods is the ease and ability to adapt 
to an evolving solution. By using Delaun- 
ay triangulation techniques, the re-gridd- 
ing in the region of interest can be 
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carried out in 0(N ) operations for 2-D 

5/2 

applications and in 0(N ) operations for 

3-D applications. 

While the unstructured-grid techniques 
offer ease of grid generation and grid 
adaptation, the solvers used in these 
methods do not appear as computationally 
efficient as the structured-grid solvers. 
Indeed, comparisons have found them slower 

by a factor of 3 to 4 for many test 
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cases (the estimate given is for a 2-D 
code here; the estimate for 3-D viscous 
flows is much worse for the same accuracy; 
see for example the timings given in Ref. 
21 for viscous calculations) . The 

relative merit of such comparisons is 
somewhat suspect, since these test 
problems have no geometrical complexity 
and thus belong naturally in the domain of 
structured-grid solvers . The unstruc- 
tured-grid methods do provide the power to 
analyze complex flow problems that are 
difficult to analyze using structured-grid 
solvers. Thus the development of unstruc- 
tured-grid technology is receiving con- 
siderable attention and is progressing 
rapidly at various laboratories. 

However, many major hurdles remain to 
be overcome before we have available a 
good three-dimensional unstructured-grid 
solver for viscous flows. The most sig- 
nificant of these hurdles is related to 
the directional sensitivity of the viscous 
flows. While the triangulation or tetra- 
hedral domain discretization does not have 
a preferred orientation, flows with boun- 
dary layers do have directional sensitiv- 
ity (i.e, boundary layers grow normal to 
the surface) . Hence, some directivity 
needs to be introduced into the grid 
generation (non-Delaunay and hence more 
time consuming) and solution algorithms 
(background grids for turbulence, genera- 
tion of many levels of grid for multi-grid 
implementations, etc.). Thus, while the 
Euler solver implementation is rather 
straightforward for the unstructured 
meshes, the implementation for high 


Reynold number viscous flows offers a sig- 
nificant challenge. At present, this com- 
plexity seems to be the stumbling block in 
extending these methods to three-dimensio- 
nal viscous flow problems. Obviously, new 
development in this area needs to occur. 

Zonal Methods 

In zonal methods, the computational 
domain is divided into sub-domains where 
grids are patched together. A dis- 
criminating feature of these techniques 
relative to multi-block or chimera 

schemes is that the sub-domains may have 
varying degrees of latitude In the 

modelled physics. Normally, the zonal 
boundaries will be two-dimensional 
surfaces and they will have to be regular. 
The zonal approach offers the ability to 
properly model the physics through 

solution approximation valid to 
particular zones. For example, a fully 
elliptic flow problem may be solved using 
a zonal method where the flow regime may 
be approximated by a parabolic system of 
equations in a large zone and by the full 
elliptic system in a smaller zone. 
Depending on the size of these domains a 

large savings in computational resources 
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may occur. Sankar, Bharadvaj and Tsung 
use a zonal approach employing a full 
Navier-Stokes solution zone embedded in an 
outer potential flow field to study an F5 
wing and an isolated helicopter rotor in 
hover. They show a savings of roughly 50% 
in computational time over full Navier- 
Stokes solutions for similar accuracy. 

From the perspective of high-lift system 
studies, zonal approaches offer pos- 
sibilities that are yet to be fully explo- 
red. Using zonal approaches, it may be 
possible to couple structured-grids with 
unstructured meshes to develop a flexible 
approach to three-dimensional problems. 
Another application of this technique may 
be in the analysis of separated flows, 
where the thin-layer approximations may 
break down and a full Navier-Stokes solu- 


tion may be required in some zones of the 
flow field. Such an approach has the pot- 
ential for generating a computationally 
efficient and accurate prediction method. 

IV. Algorithmic Issues 

Differencing schemes 

. There are several issues related to 
algorithms for high-lift studies that need 
to ke examined. Upwind schemes which per- 
form very well for supersonic flows have 
been known to perform rather poorly for 
low subsonic flows • Central difference 
schemes, which work well for subsonic 
flows do depend to a degree on carefully 
tuned artificial dissipation to stabilize 
calculations (blended second and fourth 
order dissipations, residual smoothing, 
etc.). However, for subsonic flow cal- 
culations the central differencing al- 
gorithm is probably the most well behaved. 
The accuracy of these two schemes should 
be studied on prototype problems by sys- 
tematic application. Based on the outcome 
it may turn out that one particular solu- 
tion algorithm is more suitable than the 
other for a given configurational 
analysis. An unstructured-grid algorithm, 
for example Mavriplis" 6 , which uses 
central differencing schemes may be more 
suitable for subsonic configuration studi- 
es, while other unstructured-grid solvers 
such as that by Batina 40 and Frink which 
use upwind-based schemes may be more 
suited for vortical flows. It is possible 
that no single scheme (central difference 
/ upwind scheme) will be appropriate for 
all cases. For example, consider the low 
aspect ratio, highly swept wing case. 
Here an upwind based scheme with 
controlled dissipation (e.g., TVD schemes) 
should predict the formation of vortices 
and their evolution in space more ac- 
curately than the central difference 
method where some smearing of the vor- 
ticity may occur due to added dissipation 
in most models. 


For unstructured-grid solvers, there are 
other issues related to their speed and 
accuracy that need to be fully explored, 
such as, whether vertex based or cell 
centered schemes are the most appropriate 
for the solver. The formal accuracy of 
these methods depends to a large degree on 
the particular reconstruction method chos- 
en. Cell vertex schemes are more economi- 
cal for three-dimensional (tetrahedral 
elements) unstructured grids , while cell 
centered schemes are more robust compared 
to cell vertex schemes 43 . Efficient cell- 
centered schemes for three-dimensional 
problems are possible with tetrahedral 
elements as demonstrated by Frink . There 
are also approaches which combine vertex 
based schemes with cell-centered approach 
for integration of fluxes (see Reference 
43) . The computational efficiency and 
accuracy of these approaches must also be 
examined in detail. Lomax suggests that 
there needs to be a further examination of 
special forms of structured grids in 2D 
and 3D to serve as a means for under- 
standing and evaluating unstructured grid 
solvers and their formal accuracies. 

While many of these issues will be con- 
sidered by code developers, the applied 
scientist working on high-lift system 
studies will probably be involved in 
developing methods and grids that will 
support solving flow fields around con- 
figurations with considerable geometrical 
complexities. It is quite conceivable 
that the most useful approach might be one 
that incorporates hybrid techniques. An 
example of this approach is a multi-zonal 
scheme employing hybrid computational al- 
gorithms, grid structures and/or flow 
equation models for the high-lift system 
configurations . 

Flow physics 

Progress in computational methods for 
high-lift systems strongly depends on the 
ability to model turbulence, and predict 
transition, flow separation and reattach- 


merit. For subsonic transport systems, in 
deployed high-lift situations, there may 
be cove, leading-edge slat, trailing-edge 
flap or main element flow separation. In 
some of these cases, the separated flow 
may be a massive shear layer which inter- 
acts with a boundary layer developing on 
an element downstream. The current state 
of the art in CFD does not address mas- 
sively separated flows adequately. This 
deficiency leaves CFD yielding rather im- 
potent analyses for high-lift systems as 
a result of the inability to predict where 
the flow begins to break down. 

There are also other important flow 
physics which CFD is at present unable to 
address. For example scale effects for 
high-lift systems do not show a consistent 
pattern . These anomalies are difficult 
to simulate computationally as a result of 
the significant computational resources 
required to compute flow conditions at 
flight Reynolds numbers . Relaminarization 
is a phenomenon which often occurs on the 
main element of a high-lift system in the 
influence of a deployed slat. This results 
as the flow on the main element ac- 
celerates around the leading edge due to 
extremely favorable pressure gradients at 
high Reynolds numbers . Viscous wake 
interactions is another area requiring 
further insightful studies. 

There are also areas that require im- 
mediate attention from a computational 
viewpoint. For example, the status of 
turbulence modelling for aerodynamic flows 
is rather primitive. There are two types 
of problems to be addressed here. One in- 
volves a limited understanding of transi- 
tional flows and boundary layers and the 
other is the inability to properly model 
turbulence. As a result of the inability 
to predict transition, calculations are 
often run in full "laminar" or " tur- 
bulent" options for many code comparisons. 
However, most experimental data are ob- 
tained for mixed laminar/turbulent flow 
fields. Some codes such as TLNS3D can be 


run in the tripped mode. That is, if the 
location of transition from the laminar to 
turbulent state is known, the code will be 
able to compute non-zero eddy viscosity in 
the turbulent region. Even this approach 
is inherently deficient since the initial 
evolution of turbulence (low turbulent 
Reynolds numbers, R T < 500) is not properly 
modeled by existing turbulence models. 
An improved understanding of transitional 
boundary layers and transitional zonal 
modelling is clearly needed. Narasimha 
advocates using a semi-empirical approach 
to the problem. 

The other problem relates to the actual 
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turbulence modelling itself. Menter 
evaluated the performance of four popular 
turbulence closure models for flows under 
adverse pressure gradients. The Baldwin- 
Lomax, Johnson-King, Baldwin-Barth and 
Wilcox's k-U) models were implemented in an 
incompressible Reynolds-averaged Navier- 
Stokes solver (INS code) . Menter conclud- 
ed that "the three non-equilibrium models 
gave significantly better results than the 
algebraic Baldwin-Lomax model" under stro- 
ng adverse pressure gradients. Con- 
clusions, which are not discussed here, 
were also presented relative to the per- 
formance of the three models. The authors 
feel that similar rigorous studies are 
crucial to understand the performance of 
these and other proposed turbulence models 
for multi-component airfoils and wings. 

In addition, we feel for multi-element 
airfoil and wing problems, turbulence mod- 
els based on length scales are more likely 
to fail since the choice of the ap- 
propriate length scale is difficult to 
Identify. Complex turbulence models which 
do not depend on length scales may be the 
only answer. Even after solving the 
length scale problem, experience with many 
higher moment methods has been that their 
performance may not be that attractive 
considering the additional complexities 
the^ introduce, see the discussion by Lum- 
ley 9 as well as the comparison of various 



turbulence models conducted in the AFOSR- 
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Stanford Turbulence Meeting . 

However, there are some promising new 
developments in turbulence modelling, such 
as, second-order closure'* and the^Re-Nor- 
malized Group (RNG) based models , that 
have appeared over the horizon. Results 
obtained using RNG methods are compared in 
Figure 9 {reproduced here from Reference 
52) with those using Baldwin-Lomax model 
for the RAE2822 airfoil case. Figure 9a 
presents the results obtained using RNG 
model and Figure 9b presents results ob- 
tained using Baldwin— Lomax model . Note 
that the RNG solutions improve the predic- 
tion of the shock structure relative to 
the Baldwin-Lomax model both in shock lo- 
cation and strength. This model has und- 
ergone further developments since that 
time and appears to be ready for ap- 
plication to two and three- dimensional 
solvers of multi-component configurations. 
For example, a modified RNG k- £ formula- 
tion 53 has been shown to produce excellent 
agreement for the classical backward- 
facing step problem which all of the other 
turbulence models in use have difficulty 
in predicting (see Reference 50 and the 
discussions pertaining to the backward 
facing step pp 275-283; pp 886-911) . A 
notable feature of the RNG formulation is 
that the model constants are not ad hoc 
and are derived by a consistent perturba- 
tion analysis of the Renormalised Navier- 
Stokes equations. These models may provi- 
de an avenue to improve the prediction for 
high angle of attack problems where the 
effects of turbulence are much more 
pronounced. As more and more reliance is 
placed in optimizing configurations for 
maximum lift to drag ratios or minimum 
drag, we wil be forced to examine in 
detail the agreement that these Navier- 
Stokes codes provide for integral 
quantities . 

V. Ongoing work 

The high— lift elements of NASA* s High 


Speed Research (HSR) and Advanced Sub- 
sonic Transport Program have provided the 
impetus for several significant high-lift 
efforts at Langley Research Center. The 
elements consist of a balanced experimen- 
tal and computational research program. 
Experimental work supporting the HSR 
Program involves testing a series of con- 
figurations with different types of lead- 
ing-edge high-lift devices (attached flow 
and vortex flaps) and planform variations 
(different leading edge sweeps and aspect 
ratios) . Data obtained in these tests 
include force and moment, surface pressure 
and flow visualization, A complementary 
computational program is being pursued to 
study the grid generation tools and com- 
putational methods required to analyze 
this class of vehicle. Due to the comple- 
xity of modelling the three dimensional 
high-lift system, the initial CFD effort 
concentrated on gridding and analyzing 
geometries with undeflected leading-edge 
devices. Once the cruise geometry is suc- 
cessfully analyzed the next step will be 
to analyze the high-lift configuration. 
Figure 10a is an example of an HSCT con- 
cept that was designed for a cruise Mach 
number of 3.0. This design consists of a 
blended wing body with a flattened or "pl- 
atypus" forebody. The configuration was 
analyzed by Victor Lessard of Vigyan with 
the multi-block version of CFL3D (a thin 
layer, upwind N-S code) . Figure 10b shows 
a comparison of the surface pressure dis- 
tributions obtained computationally with 
results obtained in the 8-foot Transonic 
Pressure Tunnel at NASA Langley. The com- 
parisons shown are for two cross sections 
on the configuration. The first is near 
the nose and the second is just upstream 
of the wing crank. Both comparisons show 
excellent agreement between the 
computations and experiment. The pressure 
peaks indicative of vortex flow are well 
captured by the analysis. Due to the 
agreement between theory and experiment 
obtained on this and other cruise 
geometries, we feel the method will prove 
useful for analyzing simple three- 


dimensional high-lift systems such as 
full-span attached flow flaps. 

Another CFD effort that is being pur- 
sued by Kevin Kjerstad of NASA Langley in 
support of the HSR Program includes evalu- 
ation of the three-dimensional unstructur- 
ed grid generator and Euler code developed 
by Frink, et al A ^ . The grid generation 
method is based on the "advancing front 
technique" and uses a structured back- 
ground grid to ease implementation of the 
grid generation process. An example of an 
unstructured grid generated for a generic 
high speed research configuration is shown 
in Figure 11a. The Euler solver, known as 
USM3D, is an upwind scheme developed for 
solving the three-dimensional Euler equa- 
tions on unstructured tetrahedral meshes. 
The code uses a cell-centered, finite- 
volume formulation with flux-difference 
splitting for spatial discretization. 
Experimental and USM3D results on the 
generic configuration at subsonic speeds 
for lift, drag and pitching moment are 
compared in Figure lib. Since the model 
has a sharp leading edge, the point of 
separation for the primary vortex is well 
defined and the Euler results should be 
reasonable. The comparisons in Figure 11b 
show excellent agreement between theory 
and experiment for all three quantities at 
the three angles of attack analyzed. The 
next step is to use the codes to analyze 
a three-dimensional high-lift system. The 
vortex flap concept is a good candidate 
for analysis since it has a sharp leading 
edge, hence, the separation point is known 
a priori. This configuration will be ana- 
lyzed in the near future. 

One of the efforts supporting the Ad- 
vanced Subsonic Transport Program at NASA 
Langley involves assessing the capability 
of various computational techniques for 
high-lift system application. We are cur- 
rently involved in assessing the capabili- 
ty of a structured-grid solver (TLNS3D) 
to predict the subsonic characteristics of 
a standard low-speed airfoil (GAW-1) at 


angles of attack up to C lmax . This solver 
is a transonic code with a central differ- 
encing scheme, that can be run with either 
the Baldwin-Lomax or the Johnson-King tur- 
bulence model. Experimental data for a 
GAW-1 airfoil obtained by McGhee and Beas- 
ley has been chosen for computational 
studies. This particular airfoil has geo- 
metrical characteristics (thick airfoil 
section with blunt trailing edge) which 
can pose problems for a grid sensitive 
algorithm. Tests were conducted between 
Mach numbers of 0.10 to 0.28 and angles of 
attack from -10° to 24°. The low Mach 
number data at M-0.15 has been chosen by 
us to study the robustness of the transon- 
ic structured-grid code. For the experi- 
ments, transition was fixed at 8% chord 
and the solver has the capability to com- 
pute laminar calculations up to this loca- 
tion and thereafter, switch to a turbulent 
calculation. The results presented below 
are only for the Baldwin-Lomax model since 
the performance of the Johnson-King model 
was affected by grid quality for this par- 
ticular geometry. Figure 12 shows the C p 
comparison between experimental measure- 
ments and computation using the Baldwin- 
Lomax turbulence model in TLNS3D . Com- 
putations employing the Johnson-King model 
were adversely affected by grid quality 
for this particular geometry and are not 
presented here. The agreement between the 
computed solution and experimental data is 
extremely favorable up to 8° angle of 
attack. Beginning at 12°, the computed 
pressure distributions show differences 
with experimental data at the trailing 
edge region. These differences become 
progressively worse at higher angles of 
attack. At 12° and beyond, it was noted 
by the experimenters that there was 
trailing edge flow separation which became 
progressively larger with angle of attack. 

The computed wall shear stress data as 
well as Mach contour plots (not shown 
here) do indicate trailing edge flow 
separation at 12° which becomes progres- 
sively worse at higher angles of attack. 
The degree and extent of agreement between 


experiment and theory in this separated 
zone is suspect since the included physics 
is deficient e.g, thin-layer approximation 
which breaks down in the vicinity of 
separation point (no streamwise viscous 
stress variation) and the turbulence model 
(Baldwin-Lomax model) used. The com- 
parison of sectional lift versus angle of 
attack (Figure 13) shows that the lift is 
predicted rather well by the code up to an 
angle of attack of 18°, indicating ' that 
sectional lift is insensitive to the minor 
differences in pressure distribution 
observed. The drag comparison (Figure 14) 
shows the predictions are only accurate up 
to an angle of attack of 8°. Obviously, 
minor differences that are observed in 
pressure distributions have a larger 
influence on than . Since, lift to 
drag ratio issues may dominate future 
system designs, the challenge to code 
validators is obvious. 
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Bonhaus, Anderson and Mavriplis are 
using the unstructured-grid solver of Mav- 
riplis to analyze multi-component airfoils 
for subsonic transport applications. The 
experimental data used in this comparison 
is from a Douglas four element configura- 
tion tested at the Langley LTPT tunnel. 
The computed pressure distribution over 
the elements have been compared against 
experimental data for angles of attack of 
0, 12, 18 and 20 degrees. These results 

were obtained using the Baldwin-Lomax tur- 
bulence model. The agreement between ex- 
perimental data and computations are ex- 
tremely good up to C lmax . At 0 degree ang- 
le of attack, (not shown here) , the big- 
gest difficulty is in predicting the slat 
pressure distribution in the cove region. 
At higher angles of attack, the computed 
results agree very well with measurements 
on the slat surface, while the prediction 
is off from experiment in the auxiliary 
flap as can be seen from Figure 15a. At 
20 degree angle of attack, Figure 15b, the 
disagreement is quite pronounced for both 
the main and aft flap indicating the wake 
viscous interactions are not fully cap- 


tured by the code. Figure 16 shows the 
lift vs angle of attack curve. Again, the 
predictions are in good agreement with 
experiment up to C lmax . Beyond C lmax , the 
computed lift curve shows an increase in 
lift with angle of attack, demonstrating 
the inability of the method to predict 
stall behavior. 

Previously in this paper discussions and 
examples of grids for multi-element air- 
foils have been presented. As has been 
stated, the gridding and analysis of a 
three-dimensional high-lift system is 
quite difficult. The grid generator and 
Euler code described above (Reference 41) 
were used by Dr. Mohammad Takallu of Lock- 
heed and Dr. Simha Dodbele of Vigyan to 
study a multi-element wing. The wing 
chosen was a unswept, semi-span wing 
consisting of a main element and full- 
span, double slotted flap. The surface 
grid on the configuration and part of the 
symmetry plane mesh is shown in Figure 17 . 
Even though this geometry is complex, the 
grid generation process was relatively 
straight forward. As data become avail- 
able, the results of the Euler analysis 
will be evaluated to determine the utility 
of the code for high-lift configuration 
analysis for attached flow conditions. 

Another difficulty with high-lift system 
design and analysis is the proximity of 
the ground and the effect this has on the 
flow field surrounding the configuration. 
Often there is a significant effect on 
lift due to the interference of the wing 
flow field and the ground. Figure 18 is 
an example of an unstructured grid develo- 
ped by Kyle Anderson of NASA Langley for 
a multi-component airfoil in ground ef- 
fect. The airfoil is placed in the proper 
orientation above the ground and then the 
grid is generated using Delaunay trian- 
gulation techniques. The ground is simu- 
lated by adding a zero transpiration 
boundary condition to the boundary below 
the airfoil. To analyze other ground 
heights for this same airfoils, the grid 


must then be regenerated with the airfoil 
placed in its new position. Due to the 
adaptability of unstructured grid genera- 
tion techniques, this requires relatively 
little input by the researcher. 

Ground effect analysis can also be done 
using a structured grid approach. Fig. 19 
is an example of a structured grid gen- 
erated by Dr. Steve Yaros of Langley for 
a National Aerospace Plane (NASP) type 
configuration in close proximity to the 
ground. In this case the analysis was 
done using a Navier-Stokes code; so a no 
slip boundary condition was imposed at the 
wall. Again, to analyze the configuration 
at different ground heights requires rege- 
neration of the grid. For this simple 
geometry in Figure 19 the process is stra- 
ightforward. However, for more complex 
configurations, generating new multi-block 
grids could be time consuming. 

VI . Future Plans 

Based on the discussions above it is 
possible to identify several areas where 
further coordinated work is needed. 

<i) Generation of a data base that can be 
used for developing models of flow physics 
for computer simulation. This requires 
the generation of detailed L-V data, shear 
stress, energy and fluctuation measuremen- 
ts. These measurements pose a significant 
challenge for multi-component airfoils 
where there are narrow regions of flow. 
A systematic effort should also be made 
to compare computed solutions using heuri- 
stic models with existing data bases. 

(ii) Further research is necessary to dev- 
elop three-dimensional structured and un- 
structured-grid generation techniques and 
development of hybrid (zonal) methods for 
viscous flows. Intelligent use of an ex- 
pert system may allow development of effi- 
cient three-dimensional blocked grids for 
geometrically complex configurations. 


(iii) A concerted effort should be made to 
develop new turbulence models for separat- 
ed flows and to test these models in flow 
codes. There are few and relatively poor 
performing turbulence models for separated 
flows currently available. Progress made 
in this field is bound to provide rich 
rewards. There are currently several 
good candidates such as the model of 
Wilcox 56 , RNG models, and second order 
closure models, that need to be validated. 

(iv) The need for co-ordinated efforts 
between industry and government laborator- 
ies needs to be addressed. With such vast 
areas of research to be done, a collabora- 
tive indust rial-government consortium 
would serve to reduce duplication of data 
and effort at this critical juncture. 
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(b) Slat and double slotted flap (Ref. 4) 

Figure 1: Examples of Multi— component 

Airfoil Configurations for Subsonic 
Transport Higb— lift Studies. 





Figure 2: Schematic of a low aspect 

ratio highly swept arrow wing con- 
figuration with leading and trailing 
edge segmented flaps. 


3(a) Patched grid system. 



3(b) Patched computational 
space. 
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3(d) vs a 

Figure 3: Multi-block grid system and 

comparison of theoretical and experi- 
mental lift curves for a GAW-1 multi- 
element airfoil arrangement. 



Experimental results. 
Computational predictions. 


Figure 4: vs a for a multiairfoil 

arrangement from Ref. 19. 










(b) 


Figure 8: An unstructured mesh for 

computation of a GAW-1 multi-element 
airfoil configuration in a deployed 
state. 


Figure 7 : An example of a "collar 

grid" for a cylinder intersecting a 
curved surface; (A simplified model 
of the Shuttle External Tank and 
liquid hydrogen feed line) from Ref. 
34. (a) Combined collar surface 

grid; (b) Slices of the completed 
collar grid. 
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C p Comparisons 
a *= 0 C 



Generic HSR Configuration 
Mach - 0.2 



C Comparisons 
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Figure lib: Euler comparison of 

Lift, Drag, and Pitching moment 
for a generic HSR configuration. 


Figure 12b: Comparison of C for 

P 

GAW-1 airfoil. 
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Figure 16: vs a comparisons for 

a multi— element configuration* 








Figure 18: Unstructured grid for a 

multi— component airfoil in ground 
effect study. 



Figure 19: Cross section of a struc- 

tured multi-block grid system for a 
high speed configuration in ground 
effects studies. 






